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ABSTRACT 

A  digital  computer  simulation  program  using  Digital  Simulation 
Language  is  produced  to  study  the  performance  of  a  hydrofoil  in  calm 
water.    Various  automatic  control  systems  are  studied  with  the  model 
constrained  to  the  pitch-heave-surge  mode  of  operation. 
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I.      INTRODUCTION 


A. 


COORDINATE  SYSTEMS 


Three  coordinate  systems  will  be  used  in  this  study  of  hydrofoil 
craft  dynamics.    They  are:    the  body  axes,  the  water  axes  and  the  earth 
axes.    Each  is  a  right  hand  orthogonal  system. 

The  body  axes  coordinate  system  (X,  Y,  Z)  shown  in  figure  1-1 

has  its  origin  at  the  center  of  gravity  of  the  craft  and  is  fixed  relative  to 

the  craft.    This  system  is  initially  coincident  with  the  origin  of  the  earth 
i 

axes  system,  however,  it  moves  with  the  craft  as  time  progresses.    X  is 


Poaitirt  v 


Pooltivo  R 
(Y&vrtnc  to  Stbd.) 


Positive  ? 
(nollir.5  t5  StM) 

Por'.UVa  U 


I 


,  Positive  W 

Orientation  of  the  Body  Axes  with  Respect  to  the 
Craft  and  Directions  of  Positive  Velocities 
Figure  1-1 

positive  forward,   Y  is   positive  to  starboard  and  Z  is   positive 
downward . 

The  earth  axes  coordinate  system  (X    ,  Y   ,  Z  )  is  fixed  relative  to 

£i  £j  u 

the  earths  surface.    The  origin  may  be  chosen  at  any  convenient  location, 
provided  the  body  system  origin  is  initially  at  the  same  point.    The  X^ 
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axis  lies  in  the  horizontal  plane;  it  is  initially  coincident  with,  and 


positive  in  the  same  direction  as,  the  body  X  axis.    The  Y    axis  lies  in 
the  horizontal  plane  and  is  positive  to  starboard  when  the  observer  is 
facing  the  positive  X    direction.    The  Z     axis  is  normal  to  the  horizontal 
and  is  positive  downward.    Water  surface  motion  and  crafts'  motion 
relative  to  calm  water  are  described  in  this  coordinate  system. 

The  water  axes  coordinate  system  is  aligned  with  the  relative 
velocity  vector  and  resolves  into  'lift',  'drag'  ,  and  'side  force'  directions 
All  hydrodynamic  model  test  data  is  taken  and  presented  as  plots  in  terms 
of  this  coordinate  system. 


Earth  X-Z  Vlnna 

?ODitiv«  «A 
(For  0  ••  p  m  0) 


fBody  X-Axia 


"*«. 


,mJ^ 


Pooitivo  0 
(For  •}>  «•  J/m  oj 
&»rth  X-Y  Flu* 


(Korl cental  Plana) 


Pooi 
Bocjy  Z-Axio     \ 


X 


^ivo  0  (For  0  =»  ty  m  0) 

Earth  X-Z  Piano 
(Vortical  Piano) 


VI ev  Looking  Aft 


Definition  of  Positive  Directions  for  Euler  Angles 

Figure  1-2 
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Detailed  derivations  of  the  transformation  matracies  between  coordi- 
nate systems  are  outlined  in     Ref .   1    ,  however,  only  the  results  will  be 
considered  here.    In  general,  the  body  axes  are  displaced  from  the  earth 
axes  by  the  Euler  angles,  Cf>  ,  &  ,   and  f  .    Figure  1-2  defines  the  posi- 
tive directions  for  the  Euler  angles  and  established  the  means  for  trans- 
formations between  earth  and  body  axes. 

The  hydrodynamic  force  and  moment  data  are  obtained  in  a  water 
axis  system  which  is  always  oriented  with  respect  to  the  relative  water 
velocity.    These  data  must  be  transformed  into  body  axes  so  that  their 
effects  can  be  included  in  the  equations  of  craft  motions.    In  general,  the 
water  axis  orientation  is  always  changing  with  respect  to  the  body  axes. 
This  necessitates  a  transformation  to  resolve  lift,  drag,  and  side  force 
quantities  into  body  axes.    The  orthogonal  water  axis  coordinate  system  is 
defined  in  figure  1-3.    The  angles  of  rotation  between  water  and  body  axes 
are   o<    and  (^     and  will  be  defined  later. 


LIFT 


SIDE  SLIP 


DRAG 


Figure  1-3     Definition  of  Water  Axes 
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B.         EQUATIONS  OF  MOTION 

The  development  of  the  six  equations  of  motion  for  a  hydrofoil 
follow  the  same  classical  lines  as  those  for  a  displacement  type  hull, 
This  development  is  clearly  outlined  in       Ref  1  and  4       and  therefore 
will  not  be  done  here.    The  six  equations  for  the  hydrofoil  are: 

Translation 


11  =  ™FX-QW  +  RV 


(i-i) 


\/=  vnFy  + PW-.RU 


(1-2) 


W=  mF*+Q(Jl-PV 


(1-3) 


Rotation 


P  =  ixXLL-QR(I22-Iyy)  +  (R+Qp)I 


x£ 


(1-4) 


Q--27yKRPfcri22)-(Pl-R1)ix2 


(1-5) 


R  =  4h~PQ(lyy-Ixx)-f^-P)lx; 


(1-6) 


Comparison  of  actual  and  simulated  data  show  that  the  translation 
terms  RV,  PV,  and  PW  may  be  neglected.    Similarly,  the  product  of 
inertia  quantity,  L-z,  is  less  than  10%  of  I       and  4%  less  than  I      . 
Consequently,  when  considered  in  conjunction  with  the  magnitudes  of 
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2      2 
(R+PQ) ,  (P  -R  )  and   (QR-P),   the  terms   containing  "L.     maybe  neglected 

for  this  craft. 

Given  the  equations  1-1  through  1-6  for  linear  and  rotational 

accelerations  in  body  axes ,  the  velocity  terms  in  body  axes  are  obtained 

by  merely  integrating  each  equation  one  time.    These  outputs  must  then 

be  transformed  to  earth  axes  to  be  of  useful  form.    The  transformation  of 

the  linear  equations  become: 

uE-  acos  e  c©sY*v(cosrsiNeSiN0-SiNr  S/m  0) 

■fWfcosf  Co$&Co$<p+$iN  Y  Sim  <f>)  (1-7) 


V£  -  KSin  T  Cos  e  +V(cos  r  Cos  0  +  Si ki  V Sikj  d  S//v  0) 

-fW(siN  rs/w©  C6S0-  Cos  t  S/w0) 


(1-8) 


Only  the  crafts  velocity  in  the  downward  direction  is  required  to 
compute  the  crafts   height  above  the  water  surface.    Therefore ,  the  verti- 
cal position  in  earth  axes  is  given  by: 


VX* 


WE  d-fc  (1-10) 

The  crafts   position  in  earth  axes  coordinates   is   found  from: 

xE=Juedt  VjvEdt        ,„„ 

The  Euler  angles,  (p  ,  O  ,  and  V  ,  are  also  of  interest  and  can  be 
found  by  integration  of  the  Euler  rates  0  ,  S  ,  and  *Y  .  The  rates  in 
turn  must  be  derived  from  the  body  angular  rates.    The  Euler  rates  are  not 
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easily  obtained  because  they  occur  about  axes  which  are  not  orthogonal. 
This  fact  can  be  appreciated  by  referring  to  figure  1-2  and  recalling  how 
the  angles  were  defined.       j    is  a  rotation  about  the  Z     axis,    0   is  a 
rotation  about  the  Y    axis,    and   0    is  a  rotation  about  the  X  axis.    The 
results  in  the  following  set  of  equations: 

0^P-fYS|N|0  (1-12) 

e  =  QCos  0  -  RSiN  (p  (i-i3) 

Trz(QS/M0fRCos^)cob9-f(^-p)SiK  e  d-14) 

Studies  of  this  craft  show  that  even  under  the  most  drastic  abnormal 

situations,  pitch  angles  are  not  expected  to  reach  10°  and  under  normal 

operating  conditions,  roll  angles  will  not  exceed  3°.    Therefore,   the 

following  assumptions  are  justified: 

SinQ  =  9  (1-15) 

Cos  6  =         1 

The  Euler  angles  are  obtained  by  integration  of  equations    1-12, 

13,    14. 

C.        CALCULATION  OF  VELOCITY  COMPONENTS 

The  various  foil  velocity  components  must  be  determined  at  this  . 
point  in  order  to  calculate  the  angles  of  attack  and  side  slip.    The  total 
velocity  components  of  each  foil  in  body  axes  can  be  expressed  in  terms 
of  the  craft  linear  and  angular  velocities  by  considering  the  craft  geometry 
as  represented  in  figure  1-4  and  Table  1. 
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The  resultant  foil  relative  velocity  equations   become: 
Center  foil 


UC=U  +  L2CFQ 

(1-16) 

Vc-V-L*CFP+LxcrR 

(1-17) 

WC^W~LXCFQ 

(1-18) 

Port  foil 

U?-  a+-LfpFQ-LyPFR 
VP=V-  L2PFP  +  LXPFR 
Wp  =  W-LyppP-  LKPp  Q 


Star!  oard  foil 


U£ZU+L^SP^"LVSFR 
VS  =  V-LssfrP+LXSFR 

W5  =  W+LvsFP-LXSFQ 


1-19) 
1-20) 
1-21) 

1-22) 
1-23) 
1-24) 
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Figure  1-4      Definition  of  Symbology 


Key: 


©     =     point  of  application  of  foil  forces 
o     =     point  of  application  of  strut  forces 

Note: 

ScCF   =    'Sees 

ScPF  =  "ScPS  =  ScSF  =  "ScSS  =  "ScMF 

S^SF  =  "S^SS  =  "SfPF  =  _LYPS 

LZPF  =  LZSF  =  ^MF 

T  ;ari/"1      T 

ZIPS  ZSS  are   variables  and  in  general  will  not  be   equal, 
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Mid  foil 

UM=    tt*-Lihp  Q  (1-25) 

VM=  V-  LiMf  PfLXMFR  (i-26) 

These  foil  velocities  are  in  body  axis  terms  and  can  be  used  directly 
to  calculate  the  angles  of  attack  and  angles  of  side  slip. 

The  angle  of  attack  of  the  foil  is  shown  in  figure  l-5a  and  defined 
as  the  angle  whose  tangent  is  the  total  relative  velocity  between  the  water 
and  the  foil  in  the  body  Z  direction,  divided  by  the  total  relative  velocity 
between  the  water  and  foil  in  the  body  X  direction. 

o<,=  ARCTAN  -^  (1-28) 

If  the  foil  is  not  aligned  with  the  body  X  axis,  then  the  foil  has  a  fixed 
angle  of  attack  and  the  total  angle  of  attack  becomes: 

The  side  slip  angle  is  shown  in  figure  l-5b  and  can  be  represented 
by: 


$i  =4ftCSlN-7=^ 


v/Uj;+V,WWA-    '  (1-30) 


During  actual  craft  operations ,  the  angles  of        and  will  be 

less  than  10°  permitting  use  of  the  small  angle  approximation.    Also,  the 
V  and  W  velocities  can  be  considered  very  small  in  relation  to  the  U 
velocity.-   Applying  these  simplifications,  the  equations  forc<    and  y 
become: 
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Figure  1-5 
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^=  TT"  +  o(L-wd  <2-31> 

&  r    -J^L  (1-32) 


D.  EXPANSION  OF  FORCE  AND  MOMENT  TERMS 

With  the  aid  of  figure  1-4  and  Table  1 ,  the  forces  and  moments  for 
the  craft  are  easily  written.    It  is  important  to  note,  that  when  a  numerical 
value  is  substituted  for  one  of  the  alphabetical  dimensions,  the  sign  must 
be  positive  or  negative  depending  on  the  relative  position  from  the  CG. 

One  oi  the  assumptions  made  upon  commencing  simulation  was  that 
the  struts  contributed  zero  lift  and  that  the  foils  contributed  zero  side  force. 
Therefore,  all  body  Y  motion  originates  at  the  struts  and  all  body  Z  motion 
originates  at  the  foils  . 

The  force  and  moment  equations  become: 

Fy^  F*cFfFxPF-,'FXeF  +FXMF  **FXCS  +  FK  PS  +  F<SS  t)ri  Jx  +TX  d"33) 

Fy  =  Fycs+FYPS  +  FYSS  +  ™  Jy  f1"34) 

F2  =  flCF+f£PF +rZSF  +  FZriF  *mflz  (1_35) 

L-    (F£Ly)pf    +(hh)$F~(FyLt)pS--(^L2:)sS~(FlL2)LS  (1_36) 

^/.2)pFf-(FxL^Sp^F,A2)AJF  t  (FKLz)is 

^C^I)PS    +(F*  Li)ss   +  "?"x  L*f  (1-37) 

N  =  (Fy^cs  * (Vx)ps  *  foMss  d"38* 
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All  of  these  forces  have  a  hydrodynamic  origin  except  for  the 
thrust  and  gravity  terms.    The  thrust  is  always  associated  with  the  body 
axis  and  requires  no  further  expansion.    The  gravity  terms  can  be  simpli- 
fied using  the  small  angle  approximation  for  the  Euler  angles.     (See 
summary  of  equations  A-31,  32,  33.) 

The  general  expression  for  hydrodynamic  forces  in  water  axes 
coordinates  is 

Ft'  =  ^flCi  L-L^S  (1-39) 

where  L,  D,  S  represent  lift,  drag  and  side  force  respectively.    The  lift, 
drag,  and  side  force  coefficients  vary  as  functions  of  angle  of  attack, 
angle  of  side  slip,  submergence,  velocity,  flap  and  rudder  positions. 

Ideally,  mathematical  expressions  would  have  been  developed  to 
correctly  depict  the  interrelation  of  all  the  variables  which  affect  the 
force  coefficients.    However,  hydrofoil  technology  has  not  advanced  to 
a  position  which  would  yield  such  expressions.    Some  of  the  major  prob- 
lems encountered  in  deriving  a  mathematical  expression  are: 

1)  The  occurrence  of  cavitation  and  ventilation  which  are 
completely  unpredictable. 

2)  Even  more  unpredictable  is  the  cessation  of  cavitation  and 
ventilation. 

3)  The  nonlinearity  of  the  hydrodynamic  coefficients. 

4)  Lack  of  sufficient  test  data  to  completely  describe  all  of 
the  above. 
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II.      HYDROFOIL  SIMULATION 

Given  the  equations  of  motion,  the  next  step  was  to  formulate  a 
program  which  would  simulate  the  model.    The  computer  language  chosen 
was  Digital  Simulation  Language  or  DSL.    DSL  was  chosen  over  FORTRAN 
because  of  its  adaptability  to  control  problems  and  the  ease  with  which 
the  user  can  insert  engineering  relationships  with  a  minimum  amount  of 
programming. 

This  computer  simulation  is  designed  for  the  HIGH  POINT  PC(H)-1 
which  has  a  canard  foil-strut  configuration.      Drawings  of  the  craft  and 
the  foils  are  shown  in  figures  2-1  and  2-2.    The  after  foil  is  one  unit,  but 
in  the  equations ,  it  is  considered  to  be  divided  into  three  separate  segments, 
i.e.     mid  panel  and  port/starboard  outboard  panels. 

The  following  assumptions  were  made  prior  to  beginning  simulation: 

1 .  Craft  equations  of  motion  are  valid  only  for  the  foilborne 
mode. 

2.  Weight  of  the  craft  remains  constant. 

3.  Hydrodynamic  coefficients  are  based  on  fully  wetted  surfaces, 
i.e.    no  cavitation  or  ventilation. 

4.  No  constraint  was  placed  on  crafts  heading. 

5.  Foils  are  considered  flat  surfaces  vice  dihedral/anhedral. 

6.  The  craft  has  no  automatic  control  system. 

7.  The  craft  may  be  operated  in  calm  water  or  in  a  seaway. 

Before  the  main  program  could  be  assembled,  subroutines  had  to  be 
written  to  find  the  hydrodynamic  coefficients  of  foil  lift  and  drag,  and  strut 
drag  and  side  slip.    Using  the  simulation  curves  in  Ref.  3  as  a  guide,  the 
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curves  were  stored  as  data  points  in  the  SETUP  routine.    The  subroutines 
INTERP  and  INTRP1  are  interpolation  routines  to  obtain  the  proper  coeffi- 
cient from  the  curves.    For  INTERP,  values  of  angle  of  attack  and  flap 
angle  are  used  to  enter  the  curves  and  obtain  a  value  of  foil  drag.     Angle 
of  side  slip  and  submergence  are  used  to  enter  INTRP1  to  obtain  strut  drag. 
The  curve  for  the  lift  coefficient  of  a  foil  in  the  fully  wetted  region  is   a 
straight  line  so  this  coefficient  is  found  by  merely  solving  the  equation 
for  the  straight  line.    Curves  plotted  from  the  subroutines   are   shown  in 
Appendix  B.     Listings  of  the  subroutines  are  shown  in  the  Computer  Program 
section.    The  main  program  was  then  assembled  using  figure  2-3  as  a  refer- 
ence for  data  flow.    A  complete  listing  of  the  main  program  is  shown  in 
the  Computer  Program  section. 

A .  SIMULATION  OBJECTIVES 

The  object  of  the  simulation  was  to  commence  the  runs  with  the 

craft  foilborne  and  in  a  steady  state  condition.    Once  this  condition  was 

achieved,  the  step-response  of  the  craft  would  be  studied  by  applying 

step  functions  to  excite  motion  along  and  about  each  body  axis.    That  is, 

perturbations  should  be  introduced  to  separately  excite: 

1.  Motion  along  the  X-axis 

2  .  Motion  about  the  X-axis 

3 .  Motion  along  the  Y-axis 

4.  Motion  about  the  Y-axis 

5.  Motion  along  the  Z-axis 
6  .  Motion  about  the  Z-axis 
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After  the  step  response  was  obtained ,  a  sinusoidal  sea  would  be 
inserted  to  observe  the  crafts  motion  in  a  seaway. 

B.  PROBLEMS  ENCOUNTERED 

There  were  three  major  problems  encountered  during  the  simulation 
that  are  worthy  of  discussion:     (1)  the  unknown  stability  of  the  craft, 

(2)  specification  of  units  used  in  the  derivation  of  the  equations,  and 

(3)  numerical  quantities  for  thrust  versus  speed. 

1.         Craft  Stability 

One  of  the  first  assumptions  made  was  that  the  craft  would 
be  stable  in  the  foilborne  mode  with  no  control  system.    Starting  from  this 
assumption,  the  first  consideration  became  that  of  balancing  all  of  the 
forces  and  moments.    A  speed  of  36  knots  was  chosen  to  insure  that  cavi- 
tation did  not  exist  at  any  of  the  wetted  surfaces  and  that  the  control  sur- 
faces were  zero.    The  fixed  angle  of  attack  was  assumed  zero. 

It  was  readily  apparent  that  zero  fixed  angle  of  attack  was 
an  invalid  assumption  because  the  angle  of  attack  was  to  be  the  control- 
ling parameter  in  stabilizing  the  model.    Reference  5  listed  the  fixed 
angles  of  attack  as  0  .332  forward  and  0.0558  aft.     These  values  were 
used,  but  the  model  became  unstable  at  approximately  4.9  seconds.    Only 
slight  variations  of  angle  of  attack  above  and  below  a  critical  value  caused 
the  model  to  either  raise  out  of  the  water  or  to  sink.    This  pointed  out  the 
strong  dependence  of  both  lift  and  drag  upon  the  angle  of  attack  and  the 
problem  of  balancing  the  force  and  moment  equations .    Referring  back  to 
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equations  1-33  through  1-38,  one  can  appreciate  the  complexity  involved 


in  balancing  these  equations  as  they  are  all  interrelated.    In  order  to 
simplify  the  problem  at  all,  only  the  forces  in  the  body  X  and  Z  directions 
and  the  moment  about  the  Y  axis  were  considered.    Another  simplification 
was  that  of  considering  the  entire  forward  foil  as  a  control  surface.    The 
Angle  of  Attack  program  was  then  written  to  solve  for  the  angle  of  attack 
which  minimized  the  force  in  the  body  Z  direction  and  the  moment  about 
the  body  Y  axis .    Once  the  force  and  moment  were  minimized ,  the  program 
then  checked  the  force  in  the  body  X  direction  and  adjusted  thrust  to  mini- 
mize it.    A  listing  of  this  program,  which  is  straight  FORTRAN,  is  shown 
in  the  Computer  Program  section.    The  angles  of  attack  determined  from 
this  program  were  then  inserted  into  the  main  program.    Significant  results 
were  obtained,  however,  the  model  remained  unstable .    This  indicated 
that  there  was  no  feedback  to  dampen  the  error,  so  the  model  was  function- 
ing as  an  open  loop  system.    Considering  the  output  quantities  to  be  linear 
accelerations  and  velocities,  angular  rates  and  angular  positions,  while 
the  input  quantity  is  reference  height  above  calm  water  and  X     is  the 
measured  quantity  which  is  fed  back,  the  block  diagram  can  be  represented 
as  shown  in  figure  2-3.    No  graphs  are  showi  of  this  mode  because  of  the 
instability  of  the  model. 

Stability  of  the  model  without  a  control  system  was  questioned 
at  this  point  and  it  was  determined  that  assumption  #6,  i.e.  the  craft  has 
no  automatic  control  system,  was  invalid.    Reference  5  confirmed  the  fact 
that  the  model  was  only  stable  in  three  degrees  of  freedom  without  a  control 
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system.    Additionally,  it  was  determined  that  when  the  model  was  con- 
strained to  the  pitch-heave-surge  mode  of  operation  it  was  still  only 
marginally  stable  without  a  control  system.    At  this  point,  a  control  sys- 
tem was  added  and  it  will  be  discussed  in  detail  in  Chapter  III. 

2.  Specification  of  Units 

Considerable  confusion  resulted  from  the  interpretation  of 
the  various  references  concerning  the  units  used  in  defining  Euler  angles 
and  angular  rotations »    This  problem  has  been  alleviated  here  by  the 
specific  definitions  in  the  Table  of  Symbols. 

3.  Thrust  vs  Speed 

Since  no  graphical  or  numerical  data  was  available  to  show 
the  thrust  vs  speed  relationship,  the  computer  was  utilized  to  determine 
a  value  of  thrust  for  a  given  speed.    Once  the  Angle  of  Attack  program  had 
found  the  proper  angle  of  attack  to  minimize  the  force  and  moment,  the  pro- 
gram then  checked  the  force  in  the  body  X  direction  and  adjusted  the  thrust 
to  minimize  this  force.    This  produced  an  initial  condition  for  thrust  which 
was  then  inserted  in  the  main  program.    The  main  simulation  program  incor- 
porated a  thrust  trim  routine  which  is  shown  on  the  main  listing. 

The  thrust  trim  routine  compares  the  magnitude  of  the  force 
in  the  X  direction  with  the  thrust  and  then  tests  the  difference  with  an 
arbitrarily  small  fixed  error.    If  the  difference  is  greater  than  the  error, 
then  thrust  is  set  equal  to  the  value  necessary  to  make  the  force  in  the  X 
direction  equal  to  zero.    At  the   same  time/   the  acceleration  in  the  X 
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direction  is   also  set  equal  to  zero.    If  the  difference  is   less  than 
the  error,  then  the  thrust  is  held  constant  and  the  acceleration  in  the  X 
direction  is  allowed  to  vary. 
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III.     AUTOMATIC  CONTROL  SYSTEM 

Once  it  became  apparent  that  the  model  would  remain  unstable 
without  a  control  system,  the  only  alternative  was  to  add  the  control  sys- 
tem.   As  was  previously  stated,  the  problem  will  be  greatly  simplified  by 
constraining  the  model  to  the  pitch-heave-surge  mode  of  operation.    In 
so  doing,  the  control  system  can  be  limited  to  the  forward  foil.    It  was 
assumed  that  once  the  control  system  for  the  forward  foil  was  operational, 
it  could  be  easily  extended  to  the  after  foil. 

A.  HEIGHT  SENSOR 

For  this  control  system,  a  height  sensor  was  incorporated  which 
sensed  the  change  in  height  of  the  CG  above  calm  water.    This  change 
was  utilized  as  a  positional  feedback  signal  to  control  the  flap  angle  on 
the  forward  foil.    This  was  accomplished  in  the  simulation  by  defining  the 
reference  height  of  the  CG  for  the  particular  speed  and  subtracting  the 
actual  height  to  obtain  a  difference  called  DELH.    This  difference  was 
multipled  by  a  gain  factor  to  produce  the  variable  FLAP  defined  as: 

FLAP    =    DELH*COEFH  (3-1) 

Figure  3-1  shows  the  contents  of  the  automatic  control  block  in  the  pro- 
gram listing.    Figure  3-2  is  the  block  diagram  of  this  system  showing 
submergence  height  as  the  measured  quantity  being  fed  back  to  the 
controller. 
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The  amount  of  flap  is  governed  by  the  unbalance  of  the  equations 
of  motion.    Since  the  flap  angle  is  directly  proportional  to  both  lift  and 
drag  it  makes  a  significant  contribution  to  the  balancing  of  the  equations 
of  motion.    A  parameter  study  was  made  and  determined  that  the  optimum 
value  for  COEFH  was  0.85.    Figures  3-3,  4,   5  show  the  response  of  the 
model  to  a  positive  5  deg/sec  disturbance  applied  to  Q  for  0.2  seconds. 
It  must  be  pointed  out  that  with  the  configuration  used,  the  flap  was 
capable  of  moving  instantaneously  because  an  inertia  term  had  not  yet 
been  added. 

With  the  model  stable  using  position  feedback,  the  next  refinement 
was  that  of  utilizing  acceleration,  velocity  and  position  signals  for  feed- 
back.   It  was  realized  very  quickly  that  this  was  an  unrealistic  approach. 

B.         RATE  FEEDBACK 

With  the  model  stable  using  position  feedback,  the  next  refinement 
was  that  of  utilizing  acceleration,  velocity,  and  position  signals  for  feed- 
back.   It  was  realized  very  quickly  that  this  was  an  unrealistic  approach 
in  that  two  of  the  three  feedback  signals  were  measured  at  the  CG, 
therefore,  neglecting  any  pitch  component.    The  equation  for  flap  then 
becomes: 

FLAP    =    GAIN1  DELH-GAIN3  W-GAIN4  W  (3-2) 

where  GAIN1,  GAIN3,  and  GAIN4  are  gain  constants.    Figure  3-6  shows 
the  contents  of  the  automatic  control  block  for  this  system  and  figure  3-7 
is  the  block  diagram.    The  model  became  unstable  for  all  values  except 
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zero,    of  GAIN3   and   GAIN4 .      No   graphs   were  plotted   because  of  the 
unstability  of  the  system. 

C.  RATE  FEEDBACK  WITH  ACCELEROMETER 

This  system  was  implemented  by  using  an  accelerometer  to  sense 
motion  at  the  forward  foil.    The  accelerometer  corrects  the  body  vertical 
acceleration  by  the  pitch  rate  translated  to  the  forward  foil.    This  pro- 
duces the  true  acceleration  at  the  forward  foil.    The  equation  is  given  by: 

Azc  "  V^W  (3"3) 

Figure  3-8  shows  the  contents  of  the  automatic  control  block  for  this 
system  and  figure  3-9  shows  the  block  diagram.    This  control  system  pro- 
vided  stable  operation  as  long  as  GAIN3  was  less  than  0.02  and  GAIN4 
was  less  than  0.2.    Figures  3-10,   11,   12  are  graphs  of  the  system 
response  to  the  same  disturbance  stated  previously. 

All  three  of  the  systems  listed  so  far  were  not  physically  realizable 
in  that  the  flap  was  provided  instantaneous  response  to  any  error  signal. 
This  is  not  the  case  in  a  real  system  because  the  flap  possesses  a  finite 
amount  of  inertia  which  produces  a  time  lag  in  the  response. 

D.  RATE  FEEDBACK  WITH  REAL  POLE 

The  instantaneous  response  problem  was  taken  care  of  in  this  system 
by  adding  an  inertia  term.    This  entails  adding  a  pole  to  the  root  locus 
and  placing  it  on  the  real  axis  in  the  left  half  plane.    The  addition  of  the 
pole  helped  to  stabilize  the  system,  but  also  increased  the  overshoot  and 
and  time  lag.    The  time  lag  associated  with  the  real  pole  caused  the 
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SUMMATION  OF  FORCES 
AND  MOMENTS 


A 


/ 


AUTOMATIC  CONTROL  SYSTEM 


GAIN1=0.65 

GAIN3=0.026 

GAIN4=0.225 

HREF= 11.5 

DELH=HREF-HAC 

AZC  =  (FZ/M)-LXCF*QOfTT 

FLAP=GAIN1^DELH+-GAIN3*AZC-GAIN4<  THETA 


vj 


HYDRODYNAMIC  COEFFICIENTS 


Figure  3-  8  .  Automatic  Control  Svstem  With  an 
Ac^elerometer 


44 


fl 

\ 

^   P3 

•>•>•« 

>    cy  CD 

.£3      Ph 

O- 

X| 

I 

'C        h  h 

M           CD  0) 

t\3        H-p 

■p 

•^  rH    QJ  CD 

£ 

Jh  -h  o  e 

cti 

O    O    O  o 

rH 

Ph  Ph  «; 

Ph 

/ 

\ 

X 

<"J 

EH 

EH 

•s 

IT 

Ph 

EH 

rd 

♦s 

H 

W 

P-( 

h3 

o 

n 

o 

\*^s 

X| 

o 

U 

XI 

CD 

rH 
O 

^ 

^ 

U 

4-9 

o 
o 

\       k 

/ 

\ 

C 

rt 

,c 

-P 

•H 

& 

X 

O 

cd 

r^ 

T3 

CD 

CD 

Ph 

CD 

-P 

cd 

cd 

bO 

C 

•H 

5: 

O 

^ 

CO 

e: 

f- 

cd 

CD 

u 

•P 

ba 

CD 

cd 

e: 

•H 

o 

p 

u 

CD 

M 

rH 

o 

CD 

o 

O 

rH 

O 

PQ 

< 

• 

1 

K"\ 

CD 

^ 

r» 

r- 

fri 

•H 

Ph 

CD 

u 

X| 


45 


o 

CI 


o 


o 
o 


PITCH  RATE  VS  TIME 

Rate  Feedback  With  Accelerometer 
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PITCH  ANGLE  VS  TIME 

Rate  Feedback  With  Accelerometer 
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SUBMERGENCE  VS  TIME 

Rate  Feedback  With  Acce-lerometer 
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motion  of  the  flap  to  lag  the  motion  of  the  model  and  thus  cause  the  flap 
to  continue  moving  after  the  error  had  ceased.    It  was  also  noted  that 
when  the  acceleration  and  pitch  angle  gain  constants  became  too  large, 
the  system  went  unstable.    Comparison  of  figures  3-13,   14,   15,  with 
3-10,   11,   12,  confirms  these  statements. 

The  contents  of  the  automatic  control  block  are  shown  in  figure 
3-16.    In  the  block  diagram  of  figure  3-17,  the  real  pole  term  is  incor- 
porated into  the  controller  block . 

E.         RATE  FEEDBACK  WITH  COMPENSATOR 

In  order  to  minimize  the  effect  of  the  inertia  term,  a  compensator 
had  to  be  included  in  the  system.    The  inertia  term  caused  a  time  lag  in 
the  response,  therefore,  a  lead  compensator  had  to  be  added  to  make  the 
overall  response  approximate  that  of  the  instantaneous  response.    The 
term  that  was  actually  inserted  was  a  computer  artificiality  because  it 
took  on  the  role  of  both  real  pole  and  compensator.    The  compensator  of 
the  simulation  takes  on  the  form  (P  S+1)/(P  S+l).    In  the  root  locus, 
another  pole  was  added,  however,  it  is  so  far  to  the  left  that  it  was 
ignored . 

Comparison  of  figures  3-18,   19,  20  with  figures  3-10,   11,   12 
shows  that  the  compensator  did  reduce  the  overshoot  and  time  lag  associ- 
ated with  the  real  pole.    Close  observation  shows  that  the  response  of  the 
model  with  the  lead  compensator  closely  approaches  the  response  of  the 
model  with  an  instantaneous  flap. 
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Pitch  Rate  vs  Time 
Real  Pole  Included 
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Pitch  Angle  vs  Time 
Real  Pole  Included 
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Submergence  vs  Time 
Real  Pole  Included 
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AUTOMATIC  CONTROL  SYSTEM 


GMN1=0.65 

GAIN3=0.026 

GA1N4=0.225 

HREF=11  .5 

DELH=HREF-HACT 

AZC=(FZ/M)-LXCFJ-QDDT 

FLAP1=GAIN1    DELH+GAIN3*AZC-GAIN4- THETA 
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Figure  3-16.  Automatic  Control  Svstem  Wi^tv 
Inertia  Included 


53 


• 

A 

^   &  3- 

X| 

l 

•a       P<  h 

h          CD    CD 

cd      Hf 

£tH    CD    CD 

•p 

?H  -H    O    E 

£ 

O    O    O    O 

03 
rH 

Ph 

P«4 

Ph  <s 

t 

z1 

\ 

X 

EH 

5 

•» 

EH 

P< 

W 

a 

W 

r-\ 

EH 

P-H 

O 

w 

<< 

p 

U 

v_^ 

CJ 

<D 

a 

XI 

rH 

XI 

i-i 

/ 

^ 

O 

-P 

o 
o 

\               \ 

/■ 

\ 

CD 
X| 


cd 

•H 

-p 

CD 


■F 

•H 


E 

CD 
■f 

CO 

w 

•H 

O 
,C 
CO 

E 

cti 


A! 
o 
o 

cH 
Ph 


I 
0) 

?- 

•H 

p4 


T3 
CD 

H 
C 


54 


PITCH  RATE  VS  TIME 
Compensator  included 
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PITCH  ANGLE  VS  TIME 
Compensator  Included 
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SUBMERGENCE  VS  TIME 
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The  contents  of  the  automatic  control  block  are  as  shown  in  the 
main  simulation  program  in  the  computer  program  section.    The  block  dia- 
gram for  the  system  is  shown  in  figure  3-21. 
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III.      CONCLUSIONS  AND  RECOMMENDATIONS  FOR  FURTHER  WORK 

An  invalid  beginning  assumption,  i.e.    the  craft  is  stable  without 
a  control  system,  caused  the  loss  of  a  great  deal  of  time  early  in  the 
research.    This  was  not  a  complete  loss,  however,  as  It  necessitated 
the  writing  of  the  ANGLE  OF  ATTACK  program  to  determine  the  angles  of 
attack  and  thrust  which  balanced  the  equations  of  motion. 

The  data  obtained  from  the  ANGLE  OF  ATTACK  program  illustrated 

the  tremendous  effect  slight  changes  in  angle  of  attack  make  on  the 

i 
motion  of  the  model.    Since  these  motions  are  so  slight  and  occur  in  such 

short  time  intervals,  it  would  not  be  possible  for  manual  control.    This 

necessitates  an  automatic  control  system  that  will  provide  the  necessary 

outputs  to  operate  the  model  according  to  the  given  inputs . 

A  stable  model  was  obtained  as  a  result  of  this  research,  however, 
it  is  somewhat  limited.    Using  the  present  model  and  expanding  it  to 
include  the  after  foils ,  the  model  would  be  completely  stable  in  six  degrees 
of  freedom. 

Recommendations  for  further  study  include; 

1 .  Incorporate  a  complete  automatic  control  system  which  would 
include  all  three  foils. 

2.  Produce  a  simulation  for  the  transition  period  between  hull- 
borne  and  foilborne  conditions . 

3-.        A  study  of  the  step  response  of  the  model  in  all  modes. 

4 .  A  study  of  the  stability  of  the  model  in  both  regular  and  random 

seas . 

5.  A  study  of  the  turning  characteristics. 
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APPENDIX  A 
SUMMARY  OF  COMPLETE  EQUATIONS 


Equations  for  Accelerations  in  Boch7"  Axes 

J    r 
^  =   vnrK   f  RV-QW 

V    =  v^Fy  +PW-RU 

W   ~  vn  F2    4-QU  -PV 

P  =4[L"QR(J2H-Jyy)^RfQP)jx. 


(A-1) 

(A-2) 

(A-3) 

(A-4) 
(A-5) 
(A-6) 


Velocities    in   Bochr  Axes 


U  = 


r 


J 


adt 


V  =[vd* 

W  =  f  w  Ai 
R«,rRdt 


(A-7) 

(A-B) 

(A-9) 
(A-10) 
(A-11) 
CA-12) 
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Transformation  of  Velocities  from  Body  to  Earth  Axes 

Uc-U  tosG  cos"?  +v(cos  KrsmesiN0-sm  rcoscp) 

f  W (cosy  Sin  9  Cos  (j>  f  $>m  Y  Sin  $)  (a- .13) 

Vc  =  Uc<dss  Sim^p  f  V (cosy  cos <p  -f-Si^r  sme  Sim  0) 

+  VJ(SlN  Y  Slkie  Co$>4>  -COSY  Slti  $)  (A-14) 

viE  r  -ULSiMe  -f-Veos  e  Sjn0  +WCose  c$s0  (a- 15) 

0  =    P+  H(N9  -  ••       (A_16) 

B  =    QCOS0  -  RSlSlty  ..  (A-17) 

¥  =  (q  Sin  ^  •/-  RCos  0)  Cos  e  4  (£-P)S/a/9  (a-18) 


Positions    in  Earth   Axes 


Ye  =Jv£dt 


(A-19) 

(A-20) 

(A-21) 

(A-?2) 
(A-23) 
(A-24) 
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Expansion  o^  Force  and  Moment  Equations 


Fx  -  L  FX[F  +S  FX(;S  -f  mgx  ■*-  Tx 
Fy  =  £  Fy[S    +  vngy 
F£  =  IT  F£;F   f  rnj*     ■ 

N=Z(FyL„);s 


(A-25) 
(A-26) 

U-27) 
(A-28) 

(A-29) 

(A-30) 


Expansion  of  Gravity  Terms 


™3x 
vngy 


-mo  Sin  e 

mo  Cos  9  £lN$ 
ma  Csse  £ss^ 


(A-31) 

(A-32) 
(A-33) 


Foil  Velocity  Components  in  Body  Axes 
Center  Foil- 

Vt=  V-LkfP  +  LWpR 

Wc  =  W  -  LXCFQ 


(A-34) 

(A-35) 

(A-36) 
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Port   Foil- 
Up  =U+LjlPFQ-LypfR  (A-37) 
'    Vp  =  V-  Lapp  P+ LxPF  R  (A-38) 
■    Vfp-  W+LYPF  P-LXPFQ  (A-39) 


Starboard  Foil- 

U$-U  +  LlsF  Q-LySFR  '                  (A-40) 

Vs  =  V-L25FP-fLXSFR  *  •     (a-41) 

Ws  =  W4-LySFP-LXSF9  (A-42) 


Mid   Foli- 


um a+LiAlFQ  (A-43) 

VM  =  V-LZA,pPf  LXNFR  (A-44) 

WM=W-LXMFQ  (A-45) 


Transformation   of   Water  Particle    Orbital   Velocity 
from   Earth   to   Body   Axes 

^Wl~  UEv/l  Cos  Y  £osQ  +VBtfi&iHVCose-V/BV,i  SinO  (a-46) 

Vwi  "  U£Wt'(c£stS/.v9  S/n 0-  Sim  ?  Ccs  £)+V£Wi(siN  ?  (    • 

S//V9  S/n'0  +•  Cos VCo's0)+  W£w£  Cos  e  S/w^  (a-47) 
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-C6s  VSiw0)  -f-WEWL  C^sa  Ces  <f>  (A-48) 

Each   of   these    equations   must   be    repeated    for   each    foil. 


Relative    Velocity    Compo nents 


Ur'  -  Ui-Uyji  .  (A-49) 

V^  =  W  "Vfc/t  (A-50) 

WT;  =  W;-v/w;  (A_51) 


One  set  of  equations  for  each  foil. 
Angles  of  Attack  and  Side  Slip 

WW 


o^=rt*CTAN  "ri/u.fi   +  ^/^d'xed  (A-52)* 


$=  *«CS,M  ^fT^T^T""  (A"55)* 


Total  Relative  Velocity  at  a  Particular  Foil  or  Strut 


*  One  set  of  equations  for  each  foil.   Angles  of  attack 
and  side  slip  are  calculated  in  radians. 
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Hydrodynamic   Forces    in   v/a+er   Axes 

(Each   foil   equation   and    strut   equation   must   be    repeated 
for  each    foil    and    strut,    respectively) 


Foil   Lift 


Foil   Drag 


FL=  zPy^ApCp  U-55) 


F„  =  -k<?V%CBF       ■  (A-56) 


Strut   Drag 


Strut  Side  Force 


Transformation  of  Foil  Forces  from  Water  Axes  to  Body 
Axes 

FX;F  =:~Ff)[F  Cosc<  Cos  ?tFaF^Lioo(  (a- 59) 

F7iF--FDcFSL^$  (-60) 

FZ1F  -  -FDiFSL^c<  ^o^p-FL*lF  CoS  o<  (A-61) 
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APPENDIX  B 


Fwd  Foil  Lift  Coefficient  vs  Ancle  of  Attack 
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Mid  Foil  Lift  Coefficient  vs  time 
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Aft  Foil  Lift  Coefficient  vs  Angle  of  Attack 
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Fwd  Foil  Drag  Coefficient  vs  Angle  of  Attack 
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